Simulating UNTB

Renata Diaz

The Neutral Theory of Biodiversity in R

Remember the rules?

In each time step…

  1. A 🎲 individual from the local community dies 💀

In each time step…

  1. A 🎲 individual from the local community dies 💀
    • P(💀) = \(\frac{1}{J}\)

In each time step…

  1. A 🎲 individual from the local community dies 💀
    • P(💀) = \(\frac{1}{J}\)
  2. They are replaced by a new individual through:
    • birth 🎊 or immigration 🧭

In each time step…

  1. A 🎲 individual from the local community dies 💀
    • P(💀) = \(\frac{1}{J}\)
  2. They are replaced by a new individual through:
    • birth 🎊 or immigration 🧭
    • P(🎊) + P(🧭) = 1

In each time step…

  1. A 🎲 individual from the local community dies 💀
    • P(💀) = \(\frac{1}{J}\)
  2. They are replaced by a new individual through:
    • birth 🎊 or immigration 🧭
    • P(🎊) + P(🧭) = 1
  3. Occasionally, speciation ✨ occurs

In each time step…

  1. A 🎲 individual from the local community dies 💀
    • P(💀) = \(\frac{1}{J}\)
  2. They are replaced by a new individual through:
    • birth 🎊 or immigration 🧭
    • P(🎊) + P(🧭) = 1
  3. Occasionally, speciation ✨ occurs
    • P(✨) << 1

remotes::install_github("role-model/roleR")
remotes::install_github("role-model/roleShiny")